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Abstract 



We show how the 4 He-abundance in the early Universe can 
be used to demonstrate that macroscopic samples of neutrinos in 
thermal equilibrium are indeed distributed according to Fermi- 
Dirac statistics. 



Cosmology is an excellent testing ground for theories and models of Par- 
ticle Physics Fundamental laws of Nature as well ambitious theoret- 
ical constructs have been scrutinized and confronted against cosmological 
data. Conversely, Particle Physics provides important clues to cosmological 
problems (e.g. provides candidates to dark matter, provides explanations to 
primordial baryogenesis, etc.) A paradigm of this fruitful symbiosis is the 
momentous prediction from Big Bang Nucleosynthesis (BBN) on the number 
of light neutrino species [Q] . 

A fundamental result in relativistic quantum field theory is the celebrated 
spin-statistics theorem. Particles with half-integer spin obey Fermi-Dirac 
statistics and particles with zero or integer spin obey Bose-Einstein statis- 
tics. This theorem is deeply rooted in very basic principles such as relativis- 
tic invariance, locality and microcausality and its experimental verification 
(and/or the confirmation of its consequences) is very important. Of course, 
the experimental evidence of the spin-statistics connexion for ordinary matter 
in macroscopic samples in overwhelming. Electrical (e.g. superconductivity) 
and thermal (e.g. specifie heats) properties of metals at low temperatures 
can only be explained if electrons in matter obey Fermi-Dirac statistics. Su- 
perfluidity too, is a reflection of the spin properties of matter. The planckian 
blackbody spectrum of radiation is just another phenomenological manifes- 
tation of the spin-statistics theorem and traces back to the Bose-Einstein 
character of statistical ensembles of photons. 

Neutrinos, however, are not usually found in macroscopic samples and 
held in thermal equilibrium at a given temperature in a laboratory. Thus, 
their statistical behaviour is difficult to be experimentally established. Hence, 
a direct verification of the spin-statistics connexion in the thermodynamic 
sense is still lacking in the case of neutrinos. Sure enough, they are spin 1/2 
particles, as a large amount of data coming both from nuclear reactors and 
particle accelerators do convincingly demonstrate. It would be nice then, as is 
the case for other elementary constituents, to reveal the statistical mechanics 
of neutrinos. 

Two systems in Nature do contain a large macroscopic ensemble of neu- 
trinos in thermal equilibrium. A hot supernova core is one of them. The 
other is the early Universe. In the cosmic system, one can find a period in 
the early history of the Universe where the spin-statistics connexion of neu- 
trinos has observational consequences. In the present paper we shall see that 
it is possible to decide whether or not a macroscopic collection of neutrinos 
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in thermal equilibrium at temperature T do indeed follow Fermi-Dirac statis- 
tics. BBN and in particular the primordial helium abundance Yp provides 
the basic tool for the analysis of this question as we try to explain in the rest 
of this paper. 

A fundamental postulate of Quantum Mechanics is that systems of iden- 
tical particles are described by either symmetric or antisymmetric wave func- 
tions. Systems of particles in thermal equilibrium whose states are described 
by symmetric (antisymmetric) wavefunctions are distributed according to 
Bose (Fermi) distribution functions. Now, it is a phenomenological fact that 
electrons (and other particles with half-integer spin) are described by anti- 
symmetric wave functions (and obey Fermi-Dirac statistics when assembled 
in large numbers and held in thermal equilibrium) and photons (and other 
particles with integer spin) are described by symmetric wavefunctions (and 
follow Bose-Einstein statistics). That this should be so and not the other 
way around is the content of the Pauli spin-statistics theorem, which follows 
very generally from relativistic quantum field theory as already mentioned 
above. The actual purpose of this investigation is to realize that this is obser- 
vationally the case for neutrinos. So, our intention is not, to doubt about the 
fermionic nature of neutrinos but rather, to explicitly display in a real phys- 
ical system a characteristic phenomenon associated to their nature. Namely, 
a large collection of neutrinos in thermal equilibrium at temperature T is 
Fermi-Dirac distributed. 

The helium-4 and other light element (i.e. deuterium, helium-3 and 
lithium-7) primordial abundances have been periodically evaluated and con- 
fronted with data embracing observations of ever increasing quality. The 
most recent analysis can be found in references 3 to 9. The theoretical cal- 
culation of the primordial helium abundance is by now standard and the 
basic arguments can be read in the textbooks (see e.g. ref. 10). The key 
quantity to be computed is the fraction X n of neutrons to all nucleons as the 
Universe cools from equilibrium temperatures well above neutrino freeze-out 
(T >> 1 MeV) down to the temperature where nucleosynthesis takes place 
(T ~ 0.1 MeV). X n is found from the evolution equation 

= -A(n - p)X n + X(p - n)(l - X n ) (1) 

where X(p <-> n) are the weak rates (per nucleon) that interconvert pro- 
tons and neutrons. These rates contain the microphysics (weak interaction 
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probabilities) and the macrophysics (thermalization). It is this latter com- 
ponent that we shall manipulate. 

The general strategy to test the statistical character of neutrinos in the 
cosmic sample will be to use Bose-Einstein statistics as an alternative tem- 
plet in order to check how dependent the helium abundance actually is on 
the Fermi-Dirac distribution functions (it could well be that the choice of 
statistics were irrelevant). We do not consider here other kinds of statistics 
(e.g. parastatistics) which could also lead to sound field theories. An other 
place where the choice of statistics matters -apart from the rates X(n <-» p) 
above- is the law T = T(t) that enters eq. (1) since the expansion of the Uni- 
verse depends on the fermionic/bosonic effective degrees of freedom. Here 
we choose all three neutrino species (in the rates only the electron neutrino 
enters) to obey Bose statistics. 

In order to contrast the Fermi behaviour against an alternative behaviour 
one should solve eq. (1) using neutrinos with the right (F-D) statistics and 
find also the neutron fraction in the case of neutrinos with the wrong (B-E) 
statistics, i.e. we replace everywhere [e E l KTv + 1) _1 by [e E / KTv — 1) _1 and 
replace Pauli blocking factors 1 — /fd for neutrinos wherever they appear 
in the conventional calculation (see ref. 10) by stimulated emission factors 
1 + Jbe ■ One may check at this point the self-consistency of using B-E 
statistics for neutrinos by proving that, at high temperature, detailed balance 
is satisfied, i.e. X(n — > p) — \{p — > n), and hence equilibrium can be 
maintained. Indeed, for (m n — m p )/ 'kT « 1, with T = T 7 = T v , 

/oo 
dq q 4 f(q/nT) (2) 
-oo 



where 



f( x ) = ±(l + e- x y\l-e x y 



the upper (lower) sign for negative (positive) x. 

We have implemented the above modifications in Kawano's version [11] of 
the BBN code by Wagoner [12] to obtain the relative effect on the primordial 
abundances of 4 ife, D, 3 He and 7 Li induced by the change in statistics. We 
define the quantities, 

A* b a, (BE) ^ FD) 
a,i ai(FD) 
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where a\ = Yp, a% = [D / H]p, = [ 3 He/H]p and 04 = [ 7 Li/H]p. Our 
results are shown in figure 1. 

Being the number of light neutrino species no longer a free parameter 
(but rather, fixed by LEP data to be 3) the only free parameter in our BBN 
calculation in t] w , related to the baryon density fig of the Universe through 

tt B = 0.0036/r 2 (T/2.726)Vo (4) 

where T is the microwave background temperature today and h defines the 
Hubble parameter H = 100/i km sec -1 /Mpc. Therefore, figure 1 shows the 
quantitites Aa^/a; as a function of rj w . The effect on the helium abundance 
Yp -the one we are primarily interested in since it is the one which is more 
severely constrained by observation- is negative and almost 4%. This number 
is the net result of two competing sources: a change in the weak interaction 
rates and a change in the expansion rate. The variation due to the the Hubble 
expansion is necessarily positive because trading fermionic degrees of freedom 
for bosonic degrees of freedom implies an acceleration of the expansion rate 
and, as a consequence a larger amount of A He produced. In fact, it amounts 
to +2.3% relative increase as we have explicitly checked. The dominating 
effect lies therefore in the X(n <-> p) rates which turns out to be negative. 

We are now prepared to confront our results with observation. The basic 
observational quantity is Yp and we shall use the primordial abundances of 
other light elements (D, 3 He and Li) to help us fixing a value or range of 
values for 7710- The recent detection of the Lyman a line of deuterium in a 
Quasar Absortion System (QAS) at high redshift (z=3.32) can be interpreted 
as determining the primordial deuterium abundance: |13|, 

[D/H] P = (1.9 - 2.5) x 1(T 4 (5) 

It is an order of magnitude larger than the measurement by the Hubble Space 
Telescope of the D abundance in the local interstellar medium (LISM) Jl5| 

[D/H] USM = (1.61HD x 10- 5 (6) 

Although strictly speaking there is no contradiction between eqs. (5) and 
(6), since eq. (6) has to be interpreted as a lower bound on the primordial 
D abundance 

[D/H] P > 1.5 x 1(T 5 , (7) 
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nonetheless they cannot be made to agree when extrapolating eq. (6) deep 
into the past by using galactic chemical evolution models. 

If one takes eq. (5) seriously then a very good overall agreement between 
theory and data can be achieved with a value of 7710 — 1-6 ±0.1 derived from 



the matching of BBN and the QAS deuterium abundance measurement []16 
The dependence of [D/H]p on rj w is much stronger than it is the case for 
Yp and, hence, small variations on deuterium abundance associated to the 
change in statistics as the ones shown in fig. 1, do not lead to appreciable 
modifications in the determination of 7710- Therefore, we can use the same 
7710 range in both, Fermi-Dirac and Bose-Einstein cases. 

The above determination of 7710 in turn implies a band of values for the 
helium abundance Yp. We display our result in fig. 2. In the figure we 
also collect data associated to recent analysis of observations of the helium 
abundance. The errors shown are 1 a errors and we refer to the original 
work [ refs. 3-7] for a full discussion of their nature and for an appraisal of 
the techniques of analysis involved. We also include a weighted average (full 
circle). The gray shaded column on the right is the standard BBN prediction 
and the darker column on the left is our result for neutrinos with the wrong 
statistics (the width of the bands reflects both the spread in rj W and the 
experimental uncertainty in the neutron half life: r n = 887 ± 2.0 sec ||17|| ). 
To 1 a, the spin-statistics connexion for neutrinos is confirmed (i.e. we can 
exclude the alternative templet). 

However, we should investigate the dependence of our statement on the 
chosen value of r] 1Q . As already mentioned above, other sets of data suggest 
lower deuterium abundances, which in turn mean higher values for t] w . A 
reasonable window in the intermediate segment of r/10 values should be 2.5 
< V10 < 3. In the high end of 7/10 we may take //io=6 which is comfortably 
larger than the upper limit 7710 = 5.27 derived from matching primordial 
7 Li abundance and BBN ||. Figure 3 shows the low, intermediate and high 
domains of 77 10 - The vertical dotted lines are the 1 a limits on the observed 
(average) helium abundance (the fat error bar in fig. 2). The gray and dark- 
gray boxes correspond to standard BBN results ("right" statistics) and to 
the results of our analysis (neutrinos with "wrong" statistics), respectively. 
A comment on the interpretation of this figure may serve as a concluding 
summary on the content of the paper. 

First, note that the low part of fig. 3 is just a restatement of fig. 2. It 
best exemplifies the purpose of this work, i.e. to show that BBN has the 
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potential to decide empirically whether neutrinos verify the spin-statistics 
theorem. The middle part of the figure tells us that B-E statistics (the 
alternative templet) cannot be ruled out which, of course, means that the 
F-D statistics of neutrinos cannot be demonstrated. It shows also that, to 
la, even standard BBN can be in trouble. Better data are therefore required. 
Finally, the top part of fig. 3 informs us that standard BBN does not work 
(even at the 2a level) and hence the statistics issue is irrelevant in this case. 

A last comment is in order. No matter what the final fate of the contro- 
versial deuterium data [13,14] shall be, the avalanche of new and better data 
in the coming years will certainly permit a positive discrimination of neu- 
trino statistics. As an illustration of this point, we note that a modest 20% 
reduction of the error bars in the helium data suffices to make the distance 
between boxes in fig. 3 larger then the 1-a error interval. 
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Figure Captions 



Fig. 1. The various quantities Adj/a^ as defined in the text, as a function 
of 7710. 

Fig. 2. Confronting F-D statistics (gray band) and B-E statistics (dark- 
gray band) with He-4 data from: ref. 3 (circle), ref. 4 (diamonds) ref. 5 
(triangles), ref. 6 (squares), ref. 7 (inverted triangle). The full circle is the 
average value. Here, 7710 = 1.6 ± 0.1. 

Fig. 3. Our analysis displayed for three domains of r]i . The two vertical 
dotted lines bracket the observationally allowed He-4 abundances as given 
by the average abundance shown in fig. 2. 
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Fig. 2 
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Fig. 3 
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